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We present an analysis of the quark spectral function above and below the critical temperature 
for deconfinement performed at zero and non-zero momentum in quenched lattice QCD using clover 
improved Wilson fermions in Landau gauge. It is found that the temporal quark correlation function 
in the deconfined phase near the critical temperature is well reproduced by a two-pole ansatz for the 
spectral function. This indicates that excitation modes of the quark field have small decay rates. The 
bare quark mass and momentum dependence of the spectral function is analyzed with this ansatz. 
In the chiral limit we find that the quark spectral function has two collective modes corresponding to 
the normal and plasmino excitations in the high temperature limit. Over a rather wide temperature 
range in the deconfined phase the pole mass of these modes at zero momentum, which corresponds 
to the thermal mass of the quark, is approximately proportional to temperature. With increasing 
bare quark masses the plasmino mode gradually disappears and the spectral function is dominated 
by a single pole. We also discuss quasi-particle properties of heavy quarks in the deconfined phase. 
In the confined phase, it is found that the pole ansatz for the spectral function fails completely. 

PACS numbers: ll.lO.Wx, 12.38.Aw, 12.38.Gc, 14.65.-q, 25.75.Nq 



I. INTRODUCTION 

The asymptotic freedom of QCD tells us that mat- 
ter at extremely high temperature (T) becomes a simple 
thermodynamic system composed of weakly interacting 
quarks and gluons. Thermodynamic quantities approach 
the Stefan-Bohzmann limit for massless quarks and glu- 
ons in the high temperature hmit. Deviation from this 
limit that arise from the temperature-dependent running 
gauge couphng g can be calculated perturbatively [1, 2]. 
The excitation properties of quarks and gluons in this re- 
gion can also be analyzed using perturbative techniques. 
The fact that T dominates over all other scales allows to 
adopt the hard-thermal loop (HTL) approximation [3], 
which enables us to calculate propagators of these fields 
in a gauge independent way. It is known that in leading 
order of HTL-resummed perturbation theory mass gaps 
arise in the excitation spectra of these degrees of freedom 
which are proportional to gT. They are called thermal 
masses [2]. These excitations acquire decay widths ow- 
ing to medium effects, which are proportional to g^T. At 
sufficiently high T these widths thus are parametrically 
negligible compared to the thermal masses. In addition, 
a novel excitation, the plasmino, appears in the quark 
propagator, which has a minimum in its dispersion rela- 
tion at nonzero momentum [4]. Various discussions ad- 
dressed the origin of such a peculiar mode [5-10] and its 
phenomenological consequences [11]. 

As T is lowered, the gauge coupling g grows and per- 



* Electronic address: karsch9bnl.gov 

^Electronic address: kitazawaaphys.scl.osaka-u.ac.jp 



turbation theory eventually breaks down. At the same 
time, other scales, which are not negligible compared 
to T emerge, and make the problem more complicated. 
This is, in particular, the case for the quark-gluon plasma 
(QGP) phase near the critical temperature of deconfine- 
ment, Tc, which is analyzed experimentally in heavy ion 
collision performed currently with the Relativistic Heavy 
Ion CoUider (RHIC) [12]. To analyze the highly non- 
perturbative properties of strongly interacting matter in 
the temperature range accessible to the RHIC experi- 
ments, a non-perturbative approach to QCD, such as lat- 
tice QCD Monte Carlo simulation, is needed. 

Since the leading order HTL-resummed perturbative 
calculations predict that the decay widths of quarks and 
gluons grow faster than the thermal masses as T is low- 
ered, one may naively expect that for temperatures near 
Tc the peaks in spectral functions of these degrees of free- 
dom disappear. In other words, a quasi-particle descrip- 
tion of these excitations would be inappropriate near Tc- 
There are, however, several arguments supporting the ex- 
istence of quasi-particles with the quantum numbers of 
quarks in this region. For example, quasi-particles have 
been used successfully to describe lattice QCD results on 
the equation of states and susceptibilities [13]. The ex- 
istence of quark quasi-particles near Tc is also suggested 
by results from lattice simulations on baryon number, 
electric charge and strangeness fluctuations [14, 15]. It 
is also notable that quark number scaling of the ellip- 
tic flow observed in the RHIC experiments [16] suggests 
the existence of quark degrees of freedom in the hot and 
dense matter created at RHIC. Direct studies of excita- 
tion properties of quarks based on flrst principle lattice 
calculations may therefore help to clarify the physics be- 
hind these findings. 
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After some pioneering work on thermal properties of 
the quark propagator in lattice QCD simulations [17], 
in [18] the authors of the present paper investigated the 
spectral properties of quarks at zero momentum for two 
values of temperature in the deconfined phase, T = 1.5Tc 
and 3Tc, in quenched lattice QCD in Landau gauge. In 
this work, we use a two-pole ansatz for the quark spectral 
function, p{u,p), motivated by the structure of p{u>,p) in 
the high temperature limit. This ansatz includes contri- 
butions of the normal and plasmino modes and allows to 
analyze the importance of thermal widths of these modes. 
It is found that the Euclidean correlation function for 
quarks on the lattice is well reproduced by this ansatz. 
The result indicates that excitation modes of quarks form 
sharp peaks in p{uj,p) even near Tc. In this temperature 
range close to Tc it also is found that the form of p{uj,p) 
is similar to that in the perturbative region. It clearly is 
quite different from that at T = which is given by field 
strength renormalization Z{p) and mass function M{p) 
[19]. 

At temperatures close to Tg, the appearance of scales 
which are not negligible compared to T invalidates the 
simple picture provided by the HTL approximation, 
in addition to the failure of the perturbative treat- 
ment. Prom the analysis of simpler models, composed 
of fermions and bosons, it is known that the structure 
of p{io, p) has a nontrivial dependence on the masses of 
fermions and bosons when these masses are comparable 
to T [6, 9, 10, 20]. Examples are the spectrum of fermions 
with a mass m in QED or a Yukawa model with a mass- 
less boson [6]. In these models p{uj,p) is given by the 
HTL approximation and, in particular, receives contri- 
butions from the plasmino mode whenever T/m is large 
enough. On the other hand, p{uj,p) approaches a free 
quark spectral function without the plasmino mode in 
the low temperature T/m limit. As discussed in 
[6], and summarized in Appendix A in the present pa- 
per, the numerical result, which has been obtain at the 
one loop order, shows that p{co,p) changes continuously 
as a function of T/m between these two limits. In [18], 
we analyzed the dependence of p{u!, p) on the bare quark 
mass, mo, for each T and found that the lattice result 
are in accordance with these findings. 

The main purpose of the present study is to extend the 
analysis of [18] to lower temperatures, closer to Tc, and 
to non-zero momentum. In addition to the temperature 
values analyzed before, we performed the simulations at 
three lower temperatures above and below Tc- We show 
that the results obtained in [18] do not change qualita- 
tively even at T = 1.25Tc. On the other hand, we find 
that the two-pole approximation completely fails below 
Tc, which indicates that excitations of quark fields with a 
narrow width do not exist in the confined phase. We also 
analyze the momentum dependence of p{uj,p) and show 
that the dispersion relations of the normal and plasmino 
modes behave reasonably at finite momentum. Further- 
more, we discuss the spectral properties of the light and 
charm quarks. 



An important aspect of the analysis of the quark spec- 
tral function is that these spectral functions directly re- 
flect the symmetries of the thermal system and are sensi- 
tive to their explicit or spontaneous breaking. One can, 
for example, clearly observe the effect of chiral symmetry 
breaking in the scalar channel of the quark propagator. 
We show that the behavior of scalar channel is quite dif- 
ferent between below and above Tc, while the scalar con- 
tribution to the spectral function becomes vanishingly 
small in the chiral limit above Tc, such a behavior is not 
observed below Tc- The behavior of the quark correlation 
function on exceptional configurations [21] is reported in 
detail in Appendix. B. 

This paper is organized as follows. In the next section, 
we review the general properties of the quark spectral 
function and discuss its structure for some special cases. 
The setup of the numerical simulation is summarized in 
Sec. III. We then discuss the numerical results for the 
spectral function for T > Tc in Sec. IV and Sec.V. In 
Sec. IV, we consider the bare quark mass dependence of 
the spectral function at zero momentum. This analysis 
is extended to finite momenta in Sec. V. Section VI is 
devoted to a discussion of the quark correlation function 
below Tc- We give a brief summary in Sec. VII. In Ap- 
pendix A, we review the fermion spectral function in a 
Yukawa model. In Appendix B, the behavior of the quark 
propagator on exceptional configuration is presented. 



II. GENERAL PROPERTIES OF QUARK 
PROPAGATOR AND SPECTRAL FUNCTION 

In this section, we summarize the definition and prop- 
erties of the quark propagator and the quark spectral 
function. While the content of this section may be famil- 
iar to many readers, this section also serves to introduce 
the notation used in subsequent sections. We also review 
the forms of the spectral function in some limiting cases 
which are of relevance in later sections. 



A. Definitions 

The excitation properties of quarks in a thermal 
medium can be extracted from the imaginary-time (Mat- 
subara) quark propagator which is defined as, 

S';,{T,x;Q,y) = (T,V;:(r,x)C(0,y)), (1) 

where r is the imaginary time restricted to the interval 
< r < 1/T. Here T,- denotes the time-ordering along 
the imaginary time, and ?/'|^(t, a;) is the quark operator, 
with Greek subscripts denoting the Dirac indices, and h 
and c representing colors. The thermal average (O) is de- 
fined by {O) = {l/Z)Ti[e-^"0], where the trace is taken 
over a complete set of quantum states and Z = Tr[e~''-'^]. 
In the following analysis, we use the propagator in mo- 
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mentum space, 
5^.(t,p) = t^^ / d3xd3ye^P-(-y)5;:t(r,aj;0,y), 

C J 

(2) 

where V denotes the volume of the system and Nc = 3 
is the number of colors. Equation (2) is referred to as 
the correlation function in the following. Since S*^^ (t, p) 
is a gauge covariant quantity, gauge-fixing conditions are 
required to obtain non-zero expectation values. Equa- 
tion (2) satisfies anti-periodic boundary condition along 
the imaginary time, 

S^.{t,p)=-S^4t-1/T,p). (3) 

For zero quark chemical potential, the thermal ensem- 
ble has a charge conjugation symmetry. The thermal 
average thus satisfies (O) = {COC), with C representing 
the unitary operator for charge conjugation. This sym- 
metry leads to the following identity of the correlation 
function 

S^4r,p)=[CS^{-r,-p)C-']^^ 

= -[CS^{l/T-T,-p)C-']^^, (4) 

where C is the charge conjugation matrix; CipC = Ctp'^ . 
In the last equality of Eq. (4), we used Eq. (3). 

The Fourier transform of Eq. (2) with respect to r is 
given by 

S^4t,p) =TY,e'^"^S^,{iUn,p), (5) 

n 

with the Matsubara frequencies for fermions w„ = (2n -|- 
1)ttT. The retarded propagator in real time is obtained 
by the analytic continuation of Eq. (5) as S^^,{u),p) = 

The spectral function is related to the retarded and 
Matsubara propagators as 

Pi,u{u),p) = -(l/7r)Im5|^,(cj,p) 

= -(I/tt) (5=^(0;,^) -7°5^(a;,p)t7°)^,, (6) 

and 

Sij,„{tOJn,p) = duj -. (7) 

From Eqs. (7) and (5), one obtains the relation between 
Pn.u(.<^,P) and S^^{t,p); 

roc g(l/2-TT)a;/T 

Using the charge conjugation symmetry, one can show 
that pij,v{uj-,p) obeys the following relation, 

= [Cp^(-a;,-p)C]^^. (9) 



B. Dirac structure 

Owing to parity and rotational symmetries, the Dirac 
structure of the quark spectral function at finite temper- 
ature is in general decomposed as 

= Po{'^,P){l°)ixv - Pv{0J,p){p ■ 'y)ij.u + Pa{u,p)lfj,„, 

(10) 

where p = \p\ and p = p/p. Here, 

po(u;,p) = TrD[p(u;,p)7°]/4, (11) 
Pv(w,:p) = TrD[p(a;,p)p • 7]/4, (12) 
p,{w,p)=Trn[p{w,p)]/A, (13) 

with TrD denoting the trace over Dirac indices. Using 
Eq. (9), it is shown that pq{oj,p) is an even function and 
p^j^a{uj,p) are odd functions of 

Pq{uj,p) ^ po{-uj,p), (14) 
Pv{oj,p) = -py{-uj,p), (15) 
Ps{uJ,p) = -ps{-uj,p). (16) 

In some special cases, the Dirac structure of the spec- 
tral function can also be decomposed by using projec- 
tion operators [7, 22]. Two such examples, which are 
of relevance in the subsequent sections, are the spec- 
tral function at zero momentum and that having chiral 
symmetry. For the former case, p^{u},p) vanishes and 
p{u),p = 0) can be decomposed with the projection op- 
erators L± = (1 ± 7°)/2 as 

p{w, 0) = p^(a;)L+7° + P^(a;)L-7°, (17) 

where 

p^{co) = ^TT^[p{u;,0)j°L±] 

= po{u;,0)±ps{uj,0). (18) 

We note that gamma matrices (7") in Eq. (17) come 
directly from the definition of the quark propagator; 

S{r) = {TMr)^/j{0))_= (TWXT)V't(o))/. 

If the system is chirally symmetric, the quark propa- 
gator must anti-commute with 75, and thus Ps{u),p) van- 
ishes in Eq. (10). In this case, p{oj,p) can be decomposed 
using the projection operators P±ip) = (1 ± j^p ■ 7)/2 
as 

p{u,p) = pP(a;,p)P+(p)7° + p^(a;,p)P-(p)7°, (19) 
where 

pl{u;,p) = ^TTu[p{L^,ph°P±{p)] 

= po{LJ,p)±py{uj,p). (20) 

We note that in general p±{oj) and p^{uj,p) are neither 
even nor odd functions. Instead, the charge conjugation 
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symmetry, Eq. (9), requires the following relations for 

p^{uj) and 



p±(a;,p) =p^(-a;,p). 



(21) 
(22) 



Finally, the chirally symmetric spectral function at 
zero momentum is simply written as p(uJ, 0) — po{u>, 0)7°. 
This function can be decomposed into the forms given 
by Eq. (17) as weU as Eq. (19) with p^{uj) p^(w, 0) = 
Po{^, 0). Only in this case do p±{co) and p±{co, 0) become 
even functions of w. 

The spectral functions, p±{ui) and p^{uj,p), which 
arise in the decomposition of the spectral function 
p{u),p), represent excitations having definite quantum 
numbers corresponding to each projection. Therefore, 
the excitation properties of quarks are more apparent 
in these channels rather than in Eqs. (11)-(13). More- 
over, using the definition of the spectral function, one 
can prove that they are non-negative, p±:{oj) > and 
p^{u),p) > 0. In the analyses presented in later sections, 
we therefore mainly refer to p±{oj) and p^{uj,0) instead 
of Eqs. (11)-(13). One can also show that po{uj,p) > 0, 
while the signatures of p^{LO,p) and Ps{co,p) are not de- 
termined from the definition. 

The correlation function -S'^,/(t,p) is also decomposed 
similarly to Eq. (10), 

S^,u{t, p) 

= So{T,p){'y°)^^ - S^{t,p){p ■ 7)^^ + Ss{t,p)1^^. 

(23) 

Using the charge conjugation symmetry, one can show 
that 5o(r,p) = So{l/T - T,p), S'v,s(t,p) = -5v,s(l/T'- 
T,p). For zero momentum and in the chiral limit, the 
correlation function reduces to 

5(t,0)7" = SI\t)L+ + S^{t)L_, (24) 
5(t,p)7° = SliT,p)P+{p) + S^{t,p)P.{p), (25) 



respectively. The charge conjugation symmetry requires 
that 



Sfir) = S^il/T - t), S^{t,p) = S^il/T - T,p), 



(26) 



while (r) and 5^ (t, p) are neither symmetric nor anti- 
symmetric. Only if the system is chirally symmetric, 
5'+ (t) becomes a symmetric function. 



C. Spectral functions in some special cases 

1. Free quarks 

The retarded propagator of a free quark with Dirac 
mass m is given by 

( , rt)= ^ = , A-(p;w)7° 

^freeK^,P) U) + it] - Ep ^ U) + it] + Ep' 

(28) 

where = (w + irj,p), and A±(p;to) = [Ep ± 7°p • 
7)/(2i?p) are the projection operators for the free quark 

with Ep = sjp^ + w? . The corresponding spectral 
function is given by p{uj,p) = p^^^^{w,p)A+{p;m)'y^ + 
pf[««(a;,p)A_(p;m)7° with 



P%'^%uj,p) = S{u;tEp). 



(29) 



The projection operators A±{p;m) satisfy A±(0;m) = 

L± and A±(p;0) = P±{p)- Therefore, p±{co) of the free 
quark with zero momentum reads 



pf{uj) = 6{u)Tm), 
and Pj.{co,p) with zero quark mass 

p±{u>,p) = 5{u;tp)- 



(30) 



(31) 



2. High temperature limit 

The quark propagator at asymptotically high temper- 
ature can be calculated perturbatively. The validity of 

HTL approximation allows to obtain gaiige independent 
result within this approach [3] . The quark propagator at 
leading order in perturbation is given by 



^HTL(w,P) = [f-5^HTL(w,p)] V 



(32) 



where 



Shtl(w,p) 

= — Qo - 7° + ^ 1 - -Qo - P • 7, 33 
p \p ) p \ p \p) ) 

is the quark self-energy with thermal mass = 
(l/6)52r2 and Qo = (l/2)ln(a; + \)l{x - 1) [2]. Since 
Eq. (32) is chirally symmetric, S\r^^{u:,p) and the cor- 
responding spectral function can be decomposed using 
the projections operators f±(p). The spectral functions 
p^(w,j>) then read 

p\_{u),p) =Zn(p)<5(w t -En(p)) + Zp(p)5(a; ± Ev{p)) 



+ Pcont.(±W,p), 



(34) 



5+(r) = 57(1/T-t). 



(27) 



where Pco\\\..iy^V) represents the contribution of the con- 
tinuum taking non-zero values in the space-like region. 
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p^{u!,p) has two poles in the time-hke region at E^{p) > 
and Ep{p) > 0, which are caUed normal and (anti- 
)plasmino modes, respectively [2]. For zero momen- 
tum, £^n(0) = Ep{0) = rriT and the residues satisfy 
■^n(p) = Zp{p) = 0.5. The spectral functions for zero 
momentum thus are given by, 

= ^[6{u} - niT) + S{u) + rriT)] . (35) 

3. Effect of a non-zero Dirac mass 

In the derivation of Eq. (32), it is assumed that not 
only g ^ 1 but also T dominates over all other scales, 
where the latter condition is required for the validity of 
the HTL approximation. If T is not large enough com- 
pared to other scales, the effect of the latter shows up 
which leads to modifications of the form of the quark 
propagator even if perturbation theory is still valid. An 
example for such a scale is the Dirac mass of the quark, 
TO. The effect of m has been first investigated in [6] for 
the case of QED and for a Yukawa model composed of 
a massive fermion and a massless boson. In these mod- 
els, the fermion spectral functions for zero momentum, 
P^{lu), take simple forms in the massless and infinite 
mass limits: For to/T ^ 0, T dominates over other scales 
and p±{u)) should approach Eq. (35), 

p^{cu)^^[d{u}-mT) + 6{iv + mT)]. (36) 

However, in the opposite limit, m/T oo, the Dirac 
mass dominates over T and the spectral function ap- 
proaches that of a free quark, 

p^{oj) c^SiuTm). (37) 

By comparing Eqs. (36) and (37), it is obvious that the 
number of poles in p±((^) is different in these limits. The 
analysis performed in [6] in the one-loop approximation 
showed that these two limits are nevertheless connected 
continuously; as m/T becomes smaller, the peak corre- 
sponding to the plasmino gradually manifests itself in 
p^(w). In Appendix A, the numerical results for this 
feature in the Yukawa model and details of the formal- 
ism are summarized. Also in QCD, if the temperature 
is high enough so that the one-loop approximation for 
the quark self-energy is valid, we find the same limiting 
behavior as in Eqs. (36) and (37). At intermediate val- 
ues of ni/T a similar behavior of p±(lj) as found in the 
model calculations is therfore also expected. Using lattice 
simulations, we will show in the following that the two 
limiting forms of the spectral function are observed even 
in the non-perturbative region near Tc [18] in Sec. IV. 
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1.4404 0.033 


58 
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48 


6.136 


1.6530 0.075 
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TABLE I: Simulation parameters on lattices of size x Nt . 
The last column labeled A'excp gives the number of exceptional 
configurations (see text). 

III. SIMULATION SETUP 

In this study, we analyze the quark correlation func- 
tion, Eq. (2), using lattice QCD simulations in the 
quenched approximation, where vacuum excitations of 
the quark anti-quark pairs are neglected. For the lat- 
tice fermion, we use non-perturbatively improved clover 
Wilson fermions [23, 24]. 

We use gauge field ensembles which have been gener- 
ated and used previously by the Bielefeld group to study 
screening masses and spectral functions [25, 26]. The 
simulation parameters are summarized in Table I. We 
calculate the fermion correlation function for five values 
of the temperature, three of which are above and the 
others are below Tc. The simulation for T > Tc is per- 
formed on lattices of three different volumina, x N^., 
and lattice spacing, a, in order to check the dependence 
of the numerical result on volume and lattice spacing. A 
column labeled csw in Table I gives the parameter for the 
clover coefficient. For configurations above Tc, we have 
checked that the average of the Polyakov loop on every 
configuration is closest to the Z{3) root on the real axis. 

To estimate the value of the lattice spacing at which 
our calculations for a given T/Tc have been performed 
we use Tc ~ 300MeV. This results from determinations 
of Tc in units of the square root of the string tension, 
Tel \f5 — 0.64 [27-29] and a value for the string tension, 
^Ja ~ 465 MeV, which is extracted from studies of the 
heavy quark potential in QCD with light quarks [30, 31]. 
We note that the resulting estimate for the lattice cut-off 
has to rely on a physical scale that needs to be taken 
from a physical, i.e. unquenched QCD, calculation. 

Quark propagators have been calculated after fix- 
ing each gauge field configuration to Landau gauge, 
^^J_A^ = 0. For this we used conventional and stochas- 
tic minimization algorithms with a stopping criterion, 
{l/3)tv\d^Af'\^ < a with a = IQ-^^ By comparing 
the correlation functions calculated with stopping crite- 
ria a = 10~^^, 10~^ and 10~^, we have checked that the 
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numerical result converges well at a = 10""^^. 

In the Wilson fcrmion formulation the bare quark 
mass, TO, is related to the hopping parameter k. A naive 
formula for this relation is 

where Kc denotes the critical hopping parameter corre- 
sponding to the chiral limit. For temperatures above 
Tc, we determine Kc from the k dependence of the quark 
propagator, as will be described more precisely in Sec. IV. 
The pole of the free Wilson fermion propagator at zero 
momentum, on the other hand, is given by 

rUp = — log(l + anio). (39) 

In the following, we use Eq. (39) to define the bare quark 
mass, since we found that the a dependence of the quark 
spectral function at large bare quark mass is smaller with 
the definition Eq. (39) rather than Eq. (38). The choice of 
the definitions of the bare quark mass, however, does not 
alter almost all discussions in this paper anyway, since 
our discussions never use the precise values of the bare 
quark mass. 

For temperatures T/Tc < 1.5 and values of the hop- 
ping parameters close to Kc we observe on some gauge 
field configurations an anomalous behavior of the quark 
propagator. The appearance of such exceptional con- 
figurations in calculations with light quarks in quenched 
QCD is a well-known problem in calculations with Wilson 
fcrmions [21]. We discuss the behavior of the quark cor- 
relation function on these configurations in Appendix B. 
As summarized there, the behavior of the quark propa- 
gator on these configurations is clearly unphysical, and 
it is easy to introduce a reasonable criterion to distin- 
guish them from the normal ones. The number of con- 
figurations identified to be exceptional is given in the 
last column of Table I labeled A/excp- We excluded these 
configurations from our analysis. The number of con- 
figurations analyzed, A^conf, which does not include the 
exceptional ones, is also given in Table I. One sees from 
this table that the number of the exceptional configu- 
rations tends to increase as T is lowered and as spatial 
volume is larger. In particular, we did not observe any 
exceptional configurations for T/Tc = 3. On lattices for 
64^ X 16 and T = 1.25Tc, on the other hand, almost 
40% configurations are identified to be the exceptional 
ones. As discussed in Appendix B, this large number 
is attributed to a strong correlation in the appearance 
of exceptional configurations against the gauge update. 
We thus consider that the analysis with remaining 60% 
configurations still makes sense; see Appendix B. 

The quark correlation function, Eq. (1), is the inverse 
of the fermion matrix K = Ip— m. To evaluate it numer- 
ically on the lattice, we solve the linear equation 

V'source = -^V'result, (40) 



for a given source vector V'sourco- In this study, we use a 
wall source with momentum p 

V'source(T, x) = ^6rfi ex.p{-ip ■ x), (41) 

and construct the quark propagator, Eq. (2), from the 
solution of Eq. (40), ^result = -^''"^V'source, as 

S{t,p) = ^e'P"'Vresult(T,x) 

X 

= l^J2''''"^^~''^S{T,x;0,y), (42) 
x,y 

where the Dirac and color indices are suppressed for sim- 
plicity. The point source ipsource{T,x) = (5t,o^x,o, on the 
other hand, is the simplest choice for the source term, 
which leads to a different formula for the correlation func- 
tion, 

X X 

(43) 

Translational invariance requires that the two definitions 
for S{t,p), Eqs. (42) and (43), should give the same re- 
sult. We have confirmed that this is indeed the case 
within statistical error. It is, however, foimd that the 
statistical error obtained with Eq. (42) is notably smaller 
than that with Eq. (43) while the numerical costs are al- 
most the same for both definitions. The advantage of 
the wall source becomes more prominent on lattices with 
larger spatial volume. This behavior is understood intu- 
itively: In Eq. (42), the propagators of the quark field 
starting from various points at r = are averaged; this 
does suppress fluctuations arising from the local struc- 
ture of gauge configurations. 

In the subsequent sections we limit our analyses to two 
cases; (1) zero momentum correlators with non-zero val- 
ues of the mass, nip, and (2) finite momentum correlators 
in the chiral limit and above Tc. The correlation func- 
tion for each case is decomposed as given in Eqs. (24) 
and (25), respectively. In order to reduce the statistical 
error of (r) and (t) optimally we make use of their 
periodicity in Euclidean time, Eq. (26), and define these 
correlation functions on the lattice, for example S^{t), 
as 

5f (r),att. = iKir) + S^(l/T - r)]. (44) 

IV. QUARK PROPAGATOR ABOVE Tc AT 
ZERO MOMENTUM 

In this section, we analyze the quark spectral func- 
tion above Tc for zero momentum but with finite bare 
quark mass. As discussed in Sec. II, the quark spectral 
function at zero momentum is decomposed into p±{co) as 
in Eq. (18). In the following, we consider p^{uj), since 
/0^(w) is then immediately obtained with Eq. (21). 
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A. Lattice correlation function and fitting ansatz 

In order to extract the spectral function p+ (w) from 
the lattice correlation function using Eq. (8) we assume a 
simple ansatz for the shape of p^{i^) including few fitting 
parameters. For the fitting function, we have tried four 
ansatze, two of which are single- and two-pole ones, 



(45) 
(46) 



Here ^1,21 and E12 > 0, are fitting parameters, which 
represent the residues and positions of poles, respectively. 
The pole aX ui = —E2 in Eq. (46) corresponds to the 
plasmino mode at high temperatures, while the pole at 
u! = El is the normal one. We have also used fitting 
functions that allow for a Gaussian widths, 



Zi -{uJ-Eif 
exp ■ 



Zi 



n 



■ exp 



-[^-Ei 

n 



■ exp ■ 



(47) 
-{u; + E2)^ 



^ 2 



(48) 



where Fi 2 are additional fitting parameters correspond- 
ing to the width of each peak. 

Comparing the values of x^/dof in correlated fits based 
on Eqs. (45) and (46), we found for all parameter sets 
analyzed in the present study that x^/dof in a fit based 
on Eq. (46) is more than two orders of magnitude smaller 
than fits based on the single pole ansatz, Eq. (45). The 
pole corresponding to the plasmino mode at a; — — i?2 
therefore is an intrinsic feature of the quark propagator 
above and is needed to describe the numerical results. 
A single-pole ansatz Eq. (45) is clearly ruled out. Using 
correlated fits based on Eqs. (47) and (48), the minimal 



X always is found at Fi = F2 = 0, 



the ansatze 



reduce to Eqs. (45) and (46), respectively. The extension 
to include the Gaussian widths therefore does not modify 
the fit at all. In the following analysis, we thus use the 
two-pole ansatz Eq. (46). 

Here, we note that the above result on the Gaussian 
widths is obtained in correlated fit. We checked that if 
we use uncorrelated fits, which neglect correlations be- 
tween different t's, the Gaussian ansatze can improve 
the x^/dof especially for large bare quark masses. The 
numerical result shows that for large bare quark masses 
even the single pole ansatz with non-zero width, Eq. (47), 
including three parameters can give smaller x^/dof than 
the four parameter fit based on Eq. (46). This shows 
that there exist strong correlations between different time 
slices on the lattice, which of course is expected. 

In Fig. 1, we show the numerical results for S'^(t) latt. 
on a lattice of size 64'^ x 16 at T = 3Tc for several val- 
ues of K. From the figure one sees that the shape of 
'S'^(T)iatt. approaches a single exponential function for 
small K, while it becomes flat and symmetric as k be- 
comes larger. In the vicinity of the source, i.e. at small 




x/T 



FIG. 



1: Lattice correlation function 

|3 



at T 



3T, 

for the lattice of size 64'' x 16 with various values of k. The 
solid lines represent the fitting result with the two-pole ansatz, 
Eq. (46). Note that the upper panel shows correlation func- 
tions for the heavier quarks on a logarithmic scale and also 
includes the correlation function for the lightest quark mass, 
which also is shown again in the lower panel on a linear scale. 
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o 

T3 




0.1 



FIG. 2: Bare quark mass dependence of the x^/dof obtained 
for two-pole fits at T/Tc = 1.25, 1.5 and 3 on lattices of size 
64^ X 16. 



and large r, we see deviations from this generic picture, 
which can be attributed to distortion effects arising from 
the presence of the source. In fact, by comparing the 
correlation functions on the lattices with Nr — 12 and 
16, we find that such a distortion is clearly seen only at 
the T value closest to the source. It is thus expected that 
this deviation arises only in the vicinity of the source and 
hence is negligible in the continuum limit. 
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In order to get control over distortion effects close to 
the source, we have performed fits using points r^in < 
T < N-r — Tniin witli Tmin — 2,3,4 and 5 for Nr ~ 16, 
and Tmin = 2,3 and 4 for N^- = 12. We have checked 
that the dependence of the fitting parameters on Tmin are 
small; the fit results obtained with different Tmin coincide 
within the statistical error. In the following analysis, we 
use Tmin = 4 and 3 for Nt = 16 and 12, respectively. 

The resulting correlation functions in the two-pole 
ansatz, Eq. (46), obtained from correlated fits are shown 
in Fig. 1 as solid lines. One sees that 5+ (T)iatt. is well 
reproduced by our fitting ansatz. In Fig. 2, we show the 
bare quark mass dependence of the x^/dof on lattices 
of size 64^ x 16 and T/T^ = 1.25, 1.5 and 3, where the 
critical hopping parameter in Eq. (38) will be defined 
in the next subsection. The figure shows that x^/dof is 
of order unity around nip = 0, which means that our fit- 
ting ansatz can describe the lattice correlation function 
well for light quarks. In particular, x^/dof is less than 
unity for nip/T < 0.3 with T/T^ = 1.5 and 3. For large 
bare quark masses and close to Tc, on the other hand, 
the two-pole ansatz eventually becomes worse. 

The success of two-pole ansatz for the quark corre- 
lation functions indicates that the excitation modes of 
quarks near but above Tc are good quasi-particles with 
small decay rates similar to those found in the perturba- 
tivc region. In terms of the complex pole of the propa- 
gator, the results suggest that the positions of the poles 
would be near the real axis atu) = Ei and —E2 with small 
imaginary parts. Provided that the positions of poles of 
the quark propagator are gauge independent [32], this 
also suggests that our results on the fitting parameters 
El and E2 have small gauge dependence. 



B. Pole structure 

In Fig. 3, we show the dependence of Ei, E2 and 
Z^jiZx + Z2) on the bare quark mass for TjTc — 1.25, 
1.5, and 3 obtained from calculations on lattices of size 
64"^ X 16. Error-bars have been estimated from a Jack- 
knife analysis. The bare quark mass nip is defined in 
Eq. (39) with k,; determined by the k dependence of 
p^{oj) as described below. 

The figure shows that the ratio ^2/(^1 -I- Z2) becomes 
larger with decreasing nip and eventually reaches 0.5 ir- 
respective of T. The numerical result for each T shows 
that El = E2 is satisfied within statistical errors at this 
point. Those results show two important features of the 
structure of p+(w) at this point. First, since p^{u}) be- 
comes an even function, the quark propagator is chirally 
symmetric at this point within the statistical error (see 
Sec. II B). Due to this feature, it is natural to define the 
hopping parameter satisfying Zi = Z2 = 0.5 to be the 
critical hopping parameter, Kc- The values of Kc defined 
in this way is given in the second column of Table II 







niT /T 


T/Tc 




N„ =64 Na ^ 00 


3 


0.133997(13) 


0.875(8) 0.771(20) 


1.5 


0.134999(10) 


0.906(8) 0.800(18) 


1.25 


0.135248(10) 


0.899(12) 0.803(24) 



TABLE 11: The second column shows the critical hopping pa- 
rameter Kc determined from k dependence of the fitting func- 
tions for lattices of size 64^ x 16. The values of the thermal 
mass mr analyzed on lattice with Na = 64, and that extrap- 
olated to the infinite volume limit = 00 are also presented 
in the right columns. 



^ . We have checked that these values are consistent with 
those obtained in [25, 26] from the vanishing of the isovec- 
tor axial current. The second observation is that p^{'^) 
at K = Kc has the same form as the spectral function 
in the high temperature limit, Eq. (35). We therefore 
define the thermal mass of the quark on the lattice as 
rriT = (El + E2)/2 at k — Kc- The value of tut for each 
T with Nu = 64 is given in the third column of Table II. 
One finds that the ratio rriT/T is insensitive to T in the 
range analyzed in this work, although it becomes slightly 
larger with decreasing T , which would bo in accordance 
with the expected parametric form at high temperature, 
rriT ^ gT. 

Figure 3 also shows that the relative strength of the 
plasmino pole, Z2/{Zi + Z2), decreases with increasing 
values of the bare mass, rrip. The spectral function /)+ (w) 
thus will eventually be dominated by a single-pole. This 
result agrees with the generic observations discussed in 
Sec. lie, i.e. p+(w) approaches the spectral function of 
free quarks, Eq. (37), as the bare quarks mass becomes 
larger (see also Appendix A). The quark mass depen- 
dence of the fitting parameters at large rUp thus is rea- 
sonable. We also note that Ei has a minimum at nip > 0, 
while E2 is a monotonically increasing function. This is 
in contrast to the one-loop result in the Yukawa model, 
summarized in Appendix A, where the position of the 
peak in p^{uj) at positive energy corresponding to Ei is 
a monotonically increasing function of m/T, while the 
absolute value of that at negative energy, E2, decreases 
monotonically. The rUp dependence of Ei and E2 de- 
termined from our lattice-QCD calculations therefore is 
qualitatively different from the perturbative result. The 
non-perturbative nature of the gluon field could be re- 
sponsible for this behavior. Indeed, the minimum of Ei 



^ Clearly, our definition of Kc introduced above is not unique. Pos- 
sible alternative definitions are, for example, the value of K at 
which (1) El = E2, or at which (2) the correlation function in 

the scalar channel S's(r, 0) becomes smallest. We have checked 
that the systematic error on Kc arising from these different def- 
initions is of the same order as the statistical error on Kc given 
in Table II. 
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FIG. 4: Lattice correlation functions near the chiral limit for 
T/Tc = 1.25, 1.5 and 3. The value of rUp is rUp/T ~ 0.1, 0.05 
and 0.08 for each T, respectively. 



becomes shallower with increasing temperature and the 
slope of E2 decreases, as can be seen in Fig. 3. The 
perturbative behavior thus may well be recovered in the 
perturbative high temperature limit. 

Finally, we shall briefly discuss the T dependence of 
the magnitude of the residues Zi and Z2. In Fig. 4 we 
show the correlation function S¥'{t) near the chiral limit. 



rUp/T ~ 0.08, for T/T^ = 1.25, 1.5 and 3. The figure 
shows that the magnitude of S^{t) is insensitive to T. 
This result indicates that the magnitude of both residues, 
Zi and Z2, does not have a strong T dependence for 
T/Tc > 1.25. 



C. Beyond the chiral limit 

On the lattices above Tc, one can solve Eq. (40) in and 
even beyond the chiral limit, since chiral symmetry is not 
spontaneously broken above Tc and the numerical calcu- 
lation does not suffer from the light Nambu-Goldstone 
mode. From Eq. (39), the hopping parameter for k > Kc 
corresponds to the negative Dirac mass. In Fig. 5, we 
show nip dependence of the fitting parameters near the 
chiral limit for T/Tc = 3 ^. We plot the numerical results 
only for nip > —0.2, since the convergence of the inver- 
sion routine to solve Eq. (40) based on the BiCGStab 
algorithm starts failing there. 

If the system possesses a charge conjugation symme- 
try, the sign of the Dirac mass does not affect any observ- 
ables. One can, however, show that the roles of p+(i-^) 
and p^{i^) are exchanged when the sign of the Dirac mass 



^ We have checked that correlators other than those in the vector 
and scalar channels vanish within statistical errors even for k > 



10 





1.4 




1 .3 




1.2 


CM 




LU 


1.1 








1 


llT 






0.9 




n P. 

u.o 




0.8 




0.6 




0.4 




0.2 








1 

\,E2(-m) 


/ 

VE2(m) 

□ 


\ 1 


'Zp/(Zi+Zo) '^^^ " 

Zi/(Zi+Z2) - 





-0.2 



0.2 



0.4 



±mp/T 



FIG. 5: Fitting parameters Ei^2 and ^2/(^1 + ^2) near the 
chiral limit for T/Tc — 3. The dotted lines are £1,2 as func- 
tions of —nip, and Zi/{Zi + Z2). 



is reversed; 



M/ 



(49) 



This formula is shown by the fact that po{uJ,p) and 
Ps{u},p) are even and odd, respectively, as functions of 
the bare quark mass, and Eqs. (14) and (16). In terms 
of the fitting parameters in the two-pole ansatz Eq. (46) , 
this requires that 

Ei{±m) = E2{Tm) , and Zi{±m) = Z2(t™)- (50) 

In Fig. 5, El and E2 as functions of — rrip and Zi/{Zi + 
Z2) are shown by the dotted lines. One sees from the 
figure that Eq. (50) is approximately satisfied within the 
statistical error. This result shows that our numerical 
result behaves reasonably around the chiral limit. The 
similar result is obtained on lattices for T = l.ST^. 
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FIG. 6: Bare quark mass dependence of parameters E\, E2, 
and Z2/{Zi + Z2) at T = STc for lattices of size 64^ x 16, 
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D. Volume and lattice spacing dependence 

In order to check the dependences of our results on the 
lattice spacing and finite volume, we analyzed the quark 
propagator at T/Tc = 3, 1.5 and 1.25 on lattices with 
different a and as shown in Table I. Results for Ei, 
E2 and Z2/{Zi + Z2) obtained at T/T^ = 3 for two differ- 
ent values of the lattice cut-off and two different physical 
volumina are shown in Fig. 6. Comparing the results 
obtained on lattices with different lattice cut-off, a, but 
same physical volume, i.e. 64^ x 16 and 48^ x 12, one 
sees that any possible cut-off dependence is statistically 
not significant in our analysis. On the other hand we 



FIG. 7: Extrapolation of the thermal mass of the quark to 
the infinite volume limit for T/Tc = 3, 1.5, and 1.25. 



find a clear dependence of these quantities on the spa- 
tial volume; when comparing lattices with aspect ratios 
Nc/Nr = 3 and 4 we find that the energy levels, Ei and 
E2, drop significantly near the chiral limit as the vol- 
ume is increased. For larger values of the bare mass nip 
the volume dependence of Ei becomes small. A similar 
behavior is observed for T/T^ = 1.5 and 1.25. 

The presence of a strong volume dependence of the 
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mD[GcV] ^2/(^1+^2) 


3 

1.5 


0.13114 0.13454(3) 
0.1290 0.13540(3) 


1.625(5) 0.057(2) 
1.534(6) 0.042(2) 



TABLE III: Pole mass of the charm quark mo = Ei and 
the strength of the plasmino mode Z2/{Zi + Z2) on 48^ x 12 
lattice for T/Tc = 3 and 1.5. The parameter k for the charm 
quark is those employed in [26]. 



quark propagator is not unexpected. In fact, for the 
emergence of the thermal mass hard gluons having mo- 
menta p ~ T play a crucial role [2, 3]. However, the 
lowest non-vanishing gluon momentum on the lattice is, 
Pmm/T = 2iT{Nr/N„), which still is larger than unity 
on lattices with aspect ratio N^/Nr =4. In fact, at 
present one cannot rule out that the discretization of low 
momenta can also be responsible for the success of the 
two-pole ansatz for p^{u)), since scattering processes giv- 
ing rise to the width of quasi-particles can be suppressed 
due to the discretization of momentum. An analysis of 
quark spectral functions on lattices with even larger spa- 
tial volume, or other than periodic boundary conditions, 
is needed in the future to properly control effects of small 
momenta. 

In the present study, we estimate the thermal mass 
rriT in the V ^ 00 limit by extrapolating the results ob- 
tained for two different volumina with A^,- = 16. The 
extrapolation of niT with the ansatz for the volume de- 
pendence mT{Nr/Ncr) = TnT{0)exp{cN^/N^) for each 
temperatures is shown in Fig. 7. This extrapolation sug- 
gests that finite volume effects may still be of the order 
of 15% in our current analysis of rriT/T. The value of 
mT(0) determined from this extrapolation is depicted in 
the far right columns of Table. II. 



E. Quark mass dependence 

Here we want to discuss quasi-particle properties of 
the physical quarks, i.e. up, down and charm. So far 
we have treated the bare quark mass as a free parameter 
thus. Clearly one can discuss the properties at physical 
values of the quark masses by choosing the bare mass rup 
appropriately. Such information may be exploited for the 
understanding of the QGP phase near Tc- 

In order to discuss properties of the quark propagator 
for physical quark mass values, we first show the bare 
quark mass, nip, dependences of the fitting parameters 
El, E2 and Z2/{Zi + Z2) in physical units (GeV) in 
Fig. 8. Throughout this subsection, we use lattices of 
size 48^ X 12 for T/Tc = 1.5 and 3; these sets of gauge 
configurations are exactly those used also in the analy- 
sis of charmonia at finite T in [26]. They are therefore 
most suitable for the comparison between properties of 
quarks and charmonia analyzed there. As discussed be- 
fore, the lattice spacing dependences of these results are 



small and the figure hardly changes even if we employ 
the finer lattices of size 64"^ x 16. 

Let us first investigate the quasi-particle property of 
the charm quark. For this purpose we can use the values 
of K corresponding to the charm quark evaluated in [26] , 
which are shown in Table III. In Fig. 8, the value of nip 
corresponding to these k values is shown for each T by 
a vertical line. One sees that the values of ^2/(^1 + Z2) 
on these lines, which arc shown in the far right column 
of Table HI, are small, ^2/(^1 + Z2) < 1. This means 
that the strength of the plasmino mode is weak and the 
structure of the quark spectral function is close to that 
of free quarks, Eq. (30), with a single pole at ui = Ei. 
Therefore, it should be reasonable to regard the charm 
quarks at these temperatures as free Dirac particles with 
a Dirac mass nijj = Ei . The value of m d for each T is 
shown in the fourth column of Table HI. 

The lattice simulations suggest the existence of sharp 
peaks in the spectral function of rjc and J/tp even above 
Tc up to T = 1.5 - 2Tc [26, 33]. It is interesting to 
compare the Dirac mass of the charm quark obtained 
here with the spectral functions of charmonia. In par- 
ticular, twice the Dirac mass, 2m/), gives a threshold for 
the decay process of the charmonia, provided that the 
potential between a quark and an antiquark vanishes at 
long distances. The numerical result in [26] shows that 
the energies of the peaks corresponding to rjc and J/^ 
for T/Tc = 1.5 are m^^ ~ 3.4GeV and mj/^ ~ 3.8GeV 

These values are clearly larger than 2niD ~ 3.1GeV. 

If the confinement potential is screened completely, 
m^^ and nij/.^, thus are resonance states inside the con- 
tinuum. At least at T/Tc i$ 1-5 the heavy quark free 
energy still has a complicated structure which still shows 
a linear rise over the distance range relevant for quarko- 
nium physics [34]. This needs to be taken into account 
in any further quantitative discussion. 

Here, we note that the values of k employed in [26] 
are not accurately corresponding to the physical charm 
quark: With these parameters the masses corresponding 
to rjc in the vacuum are about 3.4GeV and 4.1GeV on 
each lattice for T/Tc = 1-5 and 3, respectively. They 
are therefore slightly larger than the experimental value. 
These parameters therefore should be interpreted as a 
guide for the charm quark. In particular, the values 
of m£) in Table III are not the exact values for the 
charm quark. It should, however, be emphasized that 
the above arguments about the comparison between 2mD 
and masses of charmonia makes sense, because the same 
hopping parameter is employed in both analyses. 

Finally, we turn to a discussion on the light quarks. 
Figure 8 shows that the conditions Z2/{Z\ + Z2) — 0.5 



We note that the parameters used to determine the physical scale 
used in [26] and the present study arc slightly different. The 
masses in physical unit quoted here take this diflference into ac- 
count for comparison, i.e. we use our value for Tc to set the 
scale. 
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FIG. 8: Dependences of Ei, E2 and Z2/{Zi + Z2) on bare quark mass nip on 48"^ x 12 lattice for T/Tc = 3 and 1.5 in the 
physical unit. The value of rup corresponding to the charm quark employed in [26] are shown by the vertical lines. 



and El ~ E2 are satisfied at rrip corresponding to the 
light quark masses, say niq ~ O.OlGeV, for each tem- 
perature. This shows that the effect of a non-zero mp is 
negligible and the quasi-particle picture for hght quarks 
is close to that obtained in the high-temperature limit, 
Eq. (35), i.e. light quarks have a thermal mass ttt-t- This 
quasi-particle property of the light quarks suggests that 
the effect of the thermal mass of light quarks should 
be taken into account when one consider quark quasi- 
particles as the basic ingredients in the studies of ther- 
modynamics [13], quarkonia and the chiral transition [35] 
of the QGP phase near T^. The value of the bare mass for 
the strange quarks, ruq ~ 0.08 GeV, on the other hand, is 
in the intermediate region between the two simple limits 
for these temperatures. 



V. QUARK PROPAGATOR AT FINITE 
MOMENTUM 

In this section, we analyze the quark spectral func- 
tion at finite momentum on lattices with size 64'^ x 16 
for T/Tc = 1.5 and 3. Throughout this section we con- 
sider the quark propagator in the chiral limit by fixing 
K — Kc, where Kc is the critical hopping parameter deter- 
mined in the previous section. The correlation function 
on the lattice is calculated using a wall source, Eq. (42), 
with momentum p. The quark propagator in the chiral 
limit is decomposed into p^{uj,p) according to Eq. (20). 



Following the same approach used in Sec. IV at zero mo- 
mentum, we adopt the two-pole ansatz 

pliuj,p) = Zi5{uj - El) + Z25{uj + E2), (51) 

and determine four parameters from a correlated fit with 
Tmin — 4. The (5-functions at to = Ei and —E2 corre- 
spond to the normal and plasmino modes, respectively. 
We found that x^/dof with this ansatz is always smaller 
than 1.5 for all momenta analyzed in this study. This 
result means that the two-pole ansatz again reproduces 
the lattice correlation function well. 

In Fig. 9, we show the momentum dependence of 
the fitting parameters Ei, E2, and Z2/{Zi + Z2) for 
T/Tc =1.5 and 3. The horizontal axis represents the mo- 
mentum on the lattice p = (1/a) sin pa. The figure shows 
that for large momentum Z2/{Zi + Z2) rapidly decreases 
and El approaches the light cone. The spectral function 
at large momentum therefore approaches that of a free 
quark, consistent with the perturbative result. One also 
finds that E2 is always smaller than Ei, in contrast to 
the results in Sec. IV. This behavior qualitatively agrees 
with the behavior of poles in the high T limit [2]. One 
also observes from Fig. 9 that E2 enters the space-like re- 
gion at high momentum. While in one-loop perturbation 
theory the plasmino mode always exists in the time-like 
region, higher order corrections could give rise to such 
behavior of the plasmino mode. At least, such behavior 
does not contradict causality. 

An interesting property of the quark propagator in 
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FIG. 9: Dependences of the fitting parameters Ei and E2 
and tiie ratio ^2/(^1 + Z2) on the lattice momentum p = 
(l/a)sin(pa) for T/Tc = 1.5 and 3. 



the high temperature limit Eq. (32) is that the disper- 
sion relation of the plasmino has a minimum at finite 
momentum. In Fig. 9, one sees that the value of E2 
at lowest non-zero momentum on our lattice, Pmin = 
2TTT{Nr/N„) ~ 1.5T, is slightly larger than that at zero 
momentum, and the existence of such a minimum is sug- 
gested but not yet confirmed in the present analysis. A 
more detailed analysis at smaller momenta would clearly 
be needed, which requires calculations on lattices with a 
larger aspect ratio N^/Nr- 

The quark spectral function at high temperatures, 
Eq. (34), has a continuum /3cont.(w,p) in the space-like re- 
gion, which physically originates from the Landau damp- 
ing. At leading order, the spectral weight of /3cont.(w,p), 
f^^dujpcont.ii^TP): becomes 0.2 at most. The success of 

the two-pole fit for p^{ui,p) without the continuum there- 
fore seems inconsistent with the perturbative result. A 
possible reason for this feature is the discretization of mo- 
menta on the lattice, since the Landau damping giving 
rise to yOcont.(w,p) can be suppressed due to the miss- 
ing momenta p < T in our current analysis. Lattices 
with much larger spatial volume are required to clarify 
this problem as well as the detailed properties of the dis- 
persion relations including the minimum of the plasmino 
mode. 

So far, we discussed the spectral functions p^{uj,p), 
assuming that the scalar channel ps{T,p) vanishes. In or- 
der to check the validity of this assumption, we show 



0.02 



0.01 



-0.01 



-0.02 




FIG. 10: Scalar channel of lattice correlation function 
S's(r, 0)iatt. at K = Kc for T/Tc — 1.5 with several values 
of p. 



in Fig. 10 the momentum dependence of the corre- 
lation function in the scalar channel, <S's(t, p)iatt., for 
T/Tc = 1.5. The figure shows that the absolute values 
of S's(t, p)iatt. are smaller than 0.005 up to p/pmin — 3 
except for r-values next to the source where they suf- 
fer from lattice artifacts. These values are more than 
two orders smaller than the typical values of 5'J^(t, p)iatt., 
and thereby being negligibly small, indeed. The figure 
shows that deviations of >5's(T,p)iatt. from zero become 
statistically significant as p increases. This is a lattice 
artefact and is expected to arise from the momentum 
dependence of the mass term in the Wilson formula- 
tion; for free Wilson fermions the mass term is given 
by M{p) = mo + r(l — cos pa) with r being the Wilson 
parameter. The fact that Ss{t,p) is still small even at 
P = 3pmin — 4.5T shows that our lattice is fine enough 
so that the effect of explicit chiral symmetry breaking, 
which arises from the Wilson term, is well suppressed. 



VI. QUARK PROPAGATOR BELOW 

Next, we analyze the quark correlation function below 
Tc- In this section, we restrict the analysis to the lattices 
of size 48'^ x 16 for simplicity. In the upper two panels of 
Fig. 11, we show the correlation function at zero momen- 
tum, S'f (r)iatt., for T/Tc = 0.5 and 0.93 and for several 
values of k. The critical hopping parameter Kc for each 
T determined in [25] is 0.13566 and 0.13558, respectively. 
For comparison, S'^(T)iatt. above Tc for T/Tc = 1.25 and 
1.5 is shown in the lower two panels. 

Before starting the discussion of results obtained be- 
low Tc, we first recapitulate the qualitative behavior of 
•S^ (T)iatt. above Tc. As we have seen in Sec. IV, the 
following two qualitative features are observed above Tc'. 



(1) S'^(T)iatt. is well reproduced by the two-pole ansatz 
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FIG. 11: Lattice correlation functions S+ (T)iatt. for several 
values of T below and above Tc on 48'^ x 16 lattices. In each 
panel, (r) latt. is depicted for four values of k; two of which 
are close to Kc- In the lower two panels for T > Tc, the 
correlation functions obtained by the two-pole ansatz is also 
shown by solid lines. 



Eq. (46). The fitting results are shown by solid fines in 
the lower two panels. (2) S'^(r)iatt. approaches a sym- 
metric function, Eq. (27), which means that 
chiral symmetry of the quark propagator is recovered 
there. 

The upper two panels in Fig. 11 clearly show that the 
behavior of S'^(T)iatt. below Tc is qualitatively differ- 
ent from those above Tc- First, S'^(r)iatt. is concave 
in the log-scale plot at t/T < 0.6 for any value of k. 
Such structure can never be reproduced by the two-pole 
ansatz Eq. (46). In fact, we have checked that the fits 
with Eq. (46) always gives unacceptably large x^/dof be- 
low Tc- Moreover, this behavior of 5'^(r)iatt. cannot be 
reproduced even if we use more than three poles with 
positive residues. Our result thus indicates the violation 
of positivity of p^{uj) below Tc, which is found also in 
the Schwingcr-Dyson approach below Tc [36] . 

The failure of the two-pole ansatz for p^{uj) indicates 
the absence of quasi-particles corresponding to sharp 
peaks in p^{uj), and this result seems consistent with 
a naive picture that quark excitations are confined below 
Tc. 

In the last sections, we discussed that the gauge de- 
pendence of our result is expected to be small, due to the 
success of the two-pole approximation and the argument 
that the position of poles of propagators is gauge inde- 
pendent. This argument breaks down below Tc, since we 
no longer conclude anything about the position of poles. 
The violation of positivity of spectral functions could be 
an artifact of a specific choice of gauge fixing condition 
[6]. The calculation of -S^ ('r)iatt. with different gauge 
fixing conditions may provide us with further clues to 
understand this problem. 

From Fig. 11, one also finds that S'^(r)iatt. below Tc 
does not approach a symmetric function as k — > Kc in 
contrast to those above Tc- This means that the quark 
propagator does not become chirally symmetric even in 
the chiral limit, which is consistent with spontaneous chi- 
ral symmetry breaking below Tc- To see this behavior 
more clearly, we show in Fig. 12 the correlation func- 
tion in the scalar channel, S's(t, 0)iatt., for several values 
of K, for T = 0.5Tc and l.ST^. The figure shows that 
S's(r, 0)iatt. for T = 0.5Tc indeed stays finite in the limit 
K ^ Kc, while that for T = l.bTc is vanishingly small in 
this limit. 

Since the difference between the correlation functions 
below and above Tc is quite remarkable, we conclude that 
the thermal modification of gluon fields during the de- 
confinement transition strongly affects also the excita- 
tion properties of quarks propagating in this background 
field. Since our fit ansatz fails in the confined phase, 
however, the detailed structure of the quark propagator 
is less clear. The comparison with the quark propaga- 
tor at T = calculated in lattice simulations [19] and 
the Schwinger-Dyson equation [36] will give us insights 
to understand the present results. Such study, however, 
is beyond the scope of the present work. 
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12: Scalar channel of lattice correlation function 
0)iatt. for T/Tc = 0.5 and 1.5. 



VII. SUMMARY 

In this publication, we analyzed the dependence of the 
quark spectral function on temperature T, bare quark 
mass m, and momentum p in quenched lattice QCD 
with Landau gauge fixing. Above Tc, we found that the 
two-pole approximations for the spectral functions in the 
projected channels, and can well repro- 

duce the lattice correlation functions. Although further 
studies on the volume dependence is needed, this result 
indicates that the excitations of quarks have small de- 
cay width even near Tc. Below Tc, on the other hand, 
the two-pole ansatz fails completely and the behavior of 
quark correlation functions indicates the violation of pos- 
itivity of the spectral function. 

By analyzing the m and p dependence of these poles, 
we confirmed that above Tc two poles, the normal and 
plasmino excitations, appear in the quark propagator. 
The mass gaps of these excitation spectra should be in- 
terpreted as thermal masses, rather than Dirac mass. It 
is found that the ratio rriT jT ~ 0.8 is insensitive to T in 
the range analyzed in this study. As the bare quark mass 
is increased, the spectral function eventually changes its 
form from that having a thermal mass to the free quark 
form. We found that the bare quark mass of the charm 
quark is close to the latter, having a single mode with a 
Dirac mass. 

All analyses of the present study are based on the 
quenched approximation. Although this approximation 



includes the leading contribution in the high temperature 
limit [2] and thus is valid at sufhciently high T, the va- 
lidity of this approximation near Tc is nontrivial. For ex- 
ample, screening of giuons due to the polarization of the 
vacuum with virtual quark antiquark pairs is neglected 
in this approximation. The coupling to possible mesonic 
excitations [26, 33, 37], which may cause interesting ef- 
fects in the spectral properties of the quark [9, 10], are 
not incorporated, either. The comparison of the quark 
propagator between quenched and full lattice QCD simu- 
lations would tell us the importance of these effects near 
Tc- It also would be interesting to calculate perturba- 
tively higher order corrections to the quasi-particle prop- 
erties of quarks [20] . Such a calculation will help to clar- 
ify the origin of differences in the mass and momentum 
dependence of the quark propagator found in our lattice 
calculation in comparison to leading order perturbation 
theory. 

There are many open questions that need to be ana- 
lyzed in more detail in future calculations. A numerical 
simulation with a large spatial volume is an important 
subject among them. As discussed in the text, the in- 
fluence of momenta smaller than the temperature is not 
properly included in the present simulations; the small- 
est non-zero momentum possible on our lattices with as- 
pect ratio 4 and periodic boundary conditions is larger 
than T. Lattices allowing for momenta less than T will 
be necessary to systematically analyze the importance of 
low momentum excitations. The existence of a minimum 
in the plasmino mode could also be confirmed in such a 
study. The calculation of the quark propagator in a dif- 
ferent gauge is also important. It will allow to check di- 
rectly the gauge dependence of the present results. These 
subjects will be studied elsewhere. 



Acknowledgments 

The lattice simulations presented in this work 
have been carried out using the cluster com- 
puters ARMINIUS@Paderborn, BEN@ECT* and 
BAMQBielefeld. M. K. is supported by a Grant-in-Aid 
for Scientific Research by Monbu-Kagakusyo of Japan 
(No. 19840037). This work has been supported in 
part by contract DE-AC02-98CH10886 with the U.S. 
Department of Energy. 



APPENDIX A: QUARK SPECTRAL FUNCTION 
IN YUKAWA MODEL 

In this appendix, we review the spectral function of 
massive fermions coupled to a massless scalar boson via 
the Yukawa coupling at finite temperature T. While the 
results given in this appendix are essentially the same as 
those first discussed in [6] , we recapitulate them to make 
this paper self-contained. The details of the analysis in 
the Yukawa models at finite T are found, for example, in 
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[6] and [9, 10]. 
We start from the Lagrangian of a Yukawa model, 

C = '^p{i^-m)tp + ^df,(t)d''(l)~gy(j)ipi;, (Al) 

where tp and ((> denote the fermion and boson operators, 
m is the fermion mass, and gy represents the Yukawa 
couphng. We neglect the mass term of the scalar boson, 
since the purpose of this analysis is a study of quarks 
coupled to massless gauge bosons. It is argued in [6] that 
the qualitative result about the fermion spectral func- 
tion hardly changes even if we promote the scalar field 
in Eq. (Al) to the U{1) gauge field. In the following, we 
call the fermion field ip the quark. 

The quark self-energy in the imaginary-time formalism 
at one-loop order is given by, 

S(iw™,p) 

(A2) 

where 5o(jw„,p) = [zw„7°-_p-7-m]"i and I?o (^t'n , p) = 
[(ii/„)^ — are the Matsubara propagators for the 

free quark and the free scalar boson, respectively, with 
cij„ = (2?! -|- l)7rT and Vn — Iu-kT . After summation over 
n and analytic continuation, one obtains the self-energy 
in the real time, E(ti;,p) — Y^(iujm'p)\i^^-^u)- 

The self-energy has an ultraviolet divergence, 

which can be removed with a standard renormalization 
[6, 9, 10]. Here we simply neglect the T- independent 
part, St=o(w,p) = limT^o 2](ti),p), which includes the 
divergence. This approximation is justified if the tem- 
perature is high enough, since the T-dependent part, 
St/o(^7P) = ~ St=o('^jP), grows rapidly and 

dominates over Et^o(u;,p) as T is raised. The T- 
dependent part Et^^oC'^jP) does not suffer from any di- 
vergences and can be calculated without renormalization. 
The spectral function at one-loop order is then given by, 



TT (tj + i77)7'' — p • 7 — m — I](a;, p) 

(A3) 

In our formalism, m and T are the only dimensionfuU 
parameters and thus they control all properties of the 
system with a fixed Yukawa coupling gy. In particular, 
the dimensionless spectral function, p — Tp{ijj/T,p/T), 
is determined uniquely for a given ratio T jm. The limit 
T/m — > corresponds to low temperature, where p(uj,p) 
approaches the spectral function of free quarks, Eq. (29). 
The opposite limit, T/m — s- oo, represents the high tem- 
perature limit. If we take — s- in this limit, piuj^p) be- 
comes that calculated in the HTL approximation Eq. (34) 
with rriT = gyT/A. With a fixed nonzero gy, the 6- 
functions in p{uj,p) become peaks having a non-zero 
width of order gyT . One can, however, check numerically 
that the qualitative structure of p{uj,p) hardly changes 
with a small Yukawa coupling 5y ^ 1. 
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FIG. 13: Spectral function p^iuj) in the Yukawa model 
Eq. (Al) for several values of m/T with a Yukawa coupling 
gy = 1- 



To compare the spectral function in the Yukawa model 
with the results in Sec. IV, we limit our attention to zero 
momentum. In this case, p{uj, 0) is decomposed as in 
Eq. (17) with projection operators L±. In the following 
we also regard m/T as the external parameter, instead of 
T/m, since this is much more convenient for the compar- 
ison with the lattice result. From the above argument, 
one expects that the spectral function p^iio) in the limit 
m/T oo approaches the free quark form, Eq. (30), 

lim p^(w)~(5(cj-m), (A4) 

m/T^oo 

while in the opposite limit m/T — > 0, p^(w) should ap- 
proach Eq. (35), 

lim p^{Lo)~\{5{ijj-mT) + 5{ijj + mT)). (A5) 

We show the numerical result for p+ (w) for several val- 
ues oi m/T in Fig. 13. A fixed Yukawa coupling gy — 1 
is employed in this calculation: We have checked that 
the qualitative feature of the numerical result does not 
change with a variation of gy over a rather wide range. 
The figure shows that p+(w) for m/T = 0.01 qualita- 
tively reproduces Eq. (A5), having two peaks around 
uj = ±gT/4:. As m/T increases, the peak at negative en- 
ergy corresponding to the plasmino ceases to exist, and 
p^{iLj) is eventually dominated by a single peak at posi- 
tive energy uj ~ m. Although in the figure the width of 
the peak at positive energy, F, grows as m/T increases, 
one can check numerically and analytically that T/m van- 
ishes in the limit m/T — > cxd. The width of the peak 
therefore is negligible in this limit, and P^{ijj) reproduces 
Eq. (A4). 

The numerical result presented in Fig. 13 shows that 
the two limits given by Eqs. (A4) and (A5), respectively, 
are connected continuously at the one-loop order. It is 
also seen that the absolute value of the position of the 
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peak at positive (negative) energy is a monotonically in- 
creasing (decreasing) function of m/T. As discussed in 
the text, this feature is quaUtatively different from that 
observed on the lattice near but above Tc. 



APPENDIX B: EXCEPTIONAL 
CONFIGURATIONS 

As mentioned in Sec. Ill, we found that the quark cor- 
relation functions S{t,p) on some gauge configurations 
for T/Tc < 1.5 behave anomalously near the chiral limit 
and at zero momentum. In this appendix, we summarize 
the behavior of ^(r, p) on such exceptional configurations 
(EC). A criterion to determine the EC used in the present 
analysis is also described. 

For the moment, we consider the correlation function 
at zero momentum in the chiral limit k — Kc on the 
lattice of size 64^ x 16 for T/Tc = 1.5 (/3 = 6.872). The 
Dirac structure of S{t, 0) is decomposed as 



Ss{r) 


= TrD 


[^(^,0)]/4, 


(Bl) 


5'v,p(t) 


= TrD 


[7p^(r,0)]/4, 


(B2) 


'S'T.piy (t) 


-Tro 


[a^,5(r,0)]/4. 


(B3) 




= TrD 


[7^755(r,0)] /4, 


(B4) 


Sp{t) 


= TrD 


[75^(t,0)] /4. 


(B5) 



In Fig. 14, we show S'v,o('''), Ss{t), Sp{t), and St,12{t) 
on all 51 configurations. Among them, seven configura- 
tions are specified as the EC which are depicted by the 
bold lines. The numbers 1-7 are also assigned to these 
lines for better identification of each configuration. The 
correlation functions obtained on the other 44 configura- 
tions are shown by thin light-blue lines, which are, how- 
ever, almost degenerated in the lower three panels. From 
the figure, one clearly sees that the behavior of Ss{t), 
5p(r), and S't,12(''') on the EC is qualitatively differ- 
ent from other normal configurations which are approx- 
imately zero in these channels. As discussed in Sec. II, 
the chiral, parity, and rotational symmetries require that 
the correlation functions in these channels vanish. The 
behavior of these functions on the EC therefore is obvi- 
ously unphysical. On the other hand, Sy.q{t) tends to 
behave moderately even on the EC. 

Since the anomalous behavior of Ss{t), S'p(t), and 
S't,12(t) on the EC is quite evident, it is easy to introduce 
a criterion to identify the EC. Here, we introduce 

Ar^ l'^r(T)p, (B6) 

for each configuration where F defines different channels 
Eqs. (Bl) - (B5), and regard the configurations satisfying 

Ar > D, (B7) 

as the exceptional ones with D being the threshold to 
be determined empirically. We show Ap and At, 12 with 
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FIG. 14: Quark correlation functions in the vector, scalar, 
pseudo-scalar and tensor channels, S'v,o(t"), <S's(''"), Sp{t) and 
<S't,i2(t") (from top to bottom) on all configuration for a lattice 
of size 64^ x 16 for T/Tc — 1.5. The exceptional configurations 
are shown by the bold lines. 



Tmin = 3 on all configurations in Fig. 15. The horizontal 
axis represents the different gauge configurations which 
are ordered according to Monte Carlo time. One sees 
that Ap and At, 12 on the EC take values more than two 
orders of magnitude larger than the typical ones on the 
normal configurations. This result means that there is 
a wide range for the choice of D in Eq. (B7), and hence 
this criterion works well in practice. Our numerical result 
shows that the criterion Eq. (B7) is most successfully ap- 
plied to the pseudoscalar channel, Ap, and successively 
AT,ij with 1 < «,j < 3. Here, we notice that our for- 
mula for S{t,p) with the wall source Eq. (42), instead of 
Eq. (43), plays a crucial role for this clear separation be- 
tween the normal and exceptional configurations. In fact, 
with the correlation function calculated with Eq. (43), 
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FIG. 15: Values of Ap and At, 12 on all configurations for 
lattice of size 64^ x 16 and T/Tc = 1-5. 
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FIG. 16: Behavior of the correlation functions Snhir) on all 
configurations for a lattice of size 64^ X 16 and T/Tc = 1.5. 
The exceptional configurations are shown by the bold lines. 



the quark correlation functions have large fluctuations 
and the range of D which separates the EC with Eq. (B7) 
becomes narrower. 

Figure 15 also shows that the EC on this set of gauge 
configurations are strongly correlated. They correspond 
to subsequent configurations in Monte Carlo time, al- 
though the separation between these gauge configura- 



tions is a few times larger than the autocorrelation length 
measured in terms of the plaquette and Polyakov loop 
correlation functions [25]. This shows that the correla- 
tion among EC is significantly stronger and leads to a 
much larger autocorrelation length. A similar result is 
obtained for T/Tc — 1.25, although in this case we ob- 
served several groups of such successive EC. 

Next, in Fig. 16 we show the correlation function rep- 
resenting the propagation between right and left handed 
quarks, Vr.l = (1/2)(1 ± 75)V', 



^rl(t) = (T,^R(r)^L(0)) = TrD[i^5(r,0)] 



-{Ss{r)+SAr)). 



(B8) 



The figure shows that Srl (t) is close to zero as it should 
be even on the EC. On the other hand, the opposite chan- 
nel5LR(T) =TrD[S(r,0)(l-75)/2] = (^sM - 5p(r))/2 
behaves anomalously on the EC, as expected from the 
behavior of 6*3, p(t) in Fig. 14. Although on the gauge 
configurations for T/Tc — 1-5 we observed the anoma- 
lous behavior only on 5'lr(t) for all configurations, we 
checked that for T/Tc = 1-25 there appear anomalous 
behaviors in both Slr{t) and S'rl(t). For T/Tc = 1.25, 
however, only one of them tends to behave anomalously 
on each EC with a few exceptions. It is also found that 
the channel having the exceptional behavior tends to be 
common in a group of configurations appearing succes- 
sively. The results presented above indicate that there 
exist topological objects on the EC that cause the anoma- 
lous behavior of the quark correlation functions. To check 
this speculation, it would be interesting to directly mea- 
sure the topological charge on each configuration. 

Finally, we remark on a relation between EC ob- 
served in the quark correlation function and those in the 
hadronic channels. By measuring the pion correlation 
function constructed from the quark correlation function 
with the point source, Eq. (43), we confirmed that the ap- 
pearance of anomalous behavior in the pion correlation 
function is limited on the EC determined with the crite- 
rion Eq. (B7). We, however, also found that the behavior 
of pion propagator seems moderate on some configura- 
tions satisfying Eq. (B7). The latter behavior may be 
attributed to the form of the lattice correlation function: 
As mentioned above, the wall source, Eq. (42), for the 
quark propagator drastically reduces the statistical fluc- 
tuations compared to the point source, and hence allows 
the clear separation of EC with a criterion like Eq. (B7). 
For the same reason, the fluctuations in the pion chan- 
nel can be large when calculated with a point source and 
such fluctuations make the identification of the EC diffi- 
cult. We thus expect that if we would calculate the pion 
correlation function with a wall source operator similar 
to that used in Eq. (42) for the quarks, there would be 
a perfect correspondence in the appearance of the excep- 
tional behavior in both correlation functions. 
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